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FOREWORD 


This  barricade  investigation  was  conducted  under  Navy  Contract  N60921-89-M- 
2742  issued  by  the  Naval  Surface  Warfare  Center  (NAVSWC).  The  time  period  for 
the  study  was  April  1988  to  April  1989. 

The  study  was  directed  by  Frank  McCleskey  of  Kilkeary,  Scott  and  Associates. 
Both  Rose  Baker  of  NAVSWC  and  Lee  Wilson  of  Advanced  Technology,  Inc.,  made 
important  contributions  to  the  study.  Primarily,  Rose  Baker  made  all  of  the 
computer  runs  and  Lee  Wilson  documented  the  results. 

This  report  has  been  reviewed  by  Michael  Swisdak  who  also  served  as  the  contract 
technical  task  coordinator  at  NAVSWC. 


Approved  by: 
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Dr.  KurtF.  MUELLER 
Head,  Energetic  Materials  Division 
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INTRODUCTION 


This  report  is  divided  into  two  parts.  Part  I  contains  a  general  study  of  barricades; 
Part  II  involves  a  specific  investigation  of  the  barricade  effects  of  hilly  terrain  in. 
Taegu,  Korea. 

Barricades  have  been  in  use  for  many  years  in  areas  where  explosive  ordnance  is 
stored.  Barricades  provide  a  means  for  stopping  fragments  produced  by  inad¬ 
vertently  detonated  ammunition  stacks.  They  may  be  artificial  (man  made)  or 
natural  obstacles  such  as  hills.  This  report  covers  the  study  of  both  an  artificial  and  a 
natural  barricade. 

All  the  numerical  results  contained  in  this  report  were  produced  by  the 
FRAGHAZ  Computer  Program.1  In  order  to  provide  the  reader  with  a  basis  forjudg¬ 
ing  the  validity  of  the  numerical  data  contained  herein,  the  general  characteristics 
and  capabilities  of  the  FRAGHAZ  Computer  Program,  together  with  the  variables 
taken  into  account,  are  discussed  in  the  following  paragraphs. 


FRAGHAZ  COMPUTER  PROGRAM  DESCRIPTION 


PURPOSE 


The  FRAGHAZ  Computer  Program  was  developed  by  the  Naval  Surface  Warfare 
Center  (NAVSWC)  for  the  Department  of  Defense  Explosives  Safety  Board  (DDESB) 
to  provide  a  method  for  predicting  the  fragment  hazards  produced  by  the  inadvertent 
detonation  of  munitions.  The  computer  program  is  designed  such  that  it  can  be  easily 
modified  to  handle  a  variety  of  specific  problems  like  the  barricade  studies  contained 
in  this  report. 


HAZARD  VOLUME 


A  relatively  simple  mathematical  model  is  used  in  the  calculation  and  accumu¬ 
lation  of  fragment  hazard  statistics.  This  hazard  volume  model  is  shaped  like  a 
narrow  piece  of  pie  with  the  sharp  edge  fixed  to  the  face  of  the  fragment-producing 
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ammunition  stack.  Down  range  is  divided  into  1 0  0  -  ft  sectors  such  that  hazard 
statistics  can  be  accumulated  in  100-ft  range  increments.  The  height  of  the  hazard 
volume  is  equal  to  the  height  of  the  target,  the  target  in  this  report  being  a  standing 
man.  The  standing  man  is  represented  as  a  three-dimensional  rectangular  parallel¬ 
epiped  randomly  located  in  the  hazard  volume.  Because  the  hazard  volume  is  shaped 
like  a  piece  of  pie,  it  diverges  as  range  increases  at  an  angle  usually  taken  as  10  deg. 
This  10  deg  angle  is  commensurate  with  the  10-deg  sectors  used  with  actual  frag¬ 
ment  pickup  in  full-scale  tests  of  ammunition  stacks.  This  fragment  pickup  data  will 
be  discussed  later  under  FRAGHAZ  COMPUTER  PROGRAM  VALIDATION. 


STACK  FRAGMENTATION  CHARACTERISTICS 


The  FRAGHAZ  Computer  Program  uses  actual  fragmentation  data  from  small- 
scale  arena  tests  that  can  be  scaled  up  to  represent  large  ammunition  stacks.  The 
small-scale  tests  may  consist  of  one  or  more  pallets  of  a  specific  munition.  When  the 
individual  munitions  of  a  pallet  are  detonated  simultaneously,  or  nearly  so,  jets  are 
produced  between  adjacent  munitions.  For  the  FRAGHAZ  program,  these  jets  are 
called  interaction  areas.  These  interaction  areas  contain  the  high-density  and  high- 
velocity  fragmentation,  which  are  used  for  safety  purposes  to  calculate  down-range 
fragmentation  hazards.  The  fragment  characteristics  necessary  to  calculate  individ¬ 
ual  trajectories  for  each  fragment  recovered  in  specified  polar-  and  azimuthal-angle 
limits  are  obtained  from  the  small-scale  arena  tests.  For  a  pallet  of  Mk  82  bombs, 
there  are  260  fragments  greater  than  300  grains  that  define  the  fragmentation  char¬ 
acteristics  of  a  single  interaction  area.  Fragments  that  weigh  less  than  300  grains  do 
not  go  to  far-field  ranges  and,  even  at  the  short  ranges,  they  seldom  possess  the  kinet¬ 
ic  energy  at  impact  to  meet  the  threshold  hazard  kinetic  energy  criterion  of  DDESB 
(58  ft-lb). 


FRAGMENT  TRAJECTORIES 


For  each  fragment  recovered  in  the  small-scale  arena  test,  one  calculates  a  com¬ 
plete  trajectory,  using  a  fourth-order  Runge-Kutta  procedure.  As  such,  the  velocity, 
displacement,  and  trajectory  angle  in  three  dimensions  can  be  calculated  at  any  point 
along  the  trajectory.  In  practice  we  are  mostly  concerned  with  the  portion(s)  of  the 
trajectory  that  lie  within  the  hazard  volume.  It  is  only  then  that  the  fragment  may 
impact  the  target.  Since  we  know  the  complete  dynamics  of  the  fragment,  we  can 
calculate  hazard  statistics  such  as  hazard  density  and  hazard  probability  of  hit. 

There  are  two  types  of  trajectories  considered:  normal  and  ricochet.  The  normal 
trajectories  have  only  one  ground  impact,  while  the  ricochet  trajectories  have  two  or 
more  ground  impacts.  The  ricochet  equations  are  based  on  tests  conducted  by  the 
Ballistic  Research  Laboratories  in  the  late  1960s.2 
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Wind  can  be  included  in  the  trajectory  calculations.  For  a  tailwind,  the  added 
range  due  to  wind  is  approximately  equal  to  the  time  of  flight  times  the  windspeed. 
Because  time  of  flight  varies  up  to  about  30  sec,  a  tailwind  of  60  ft/sec  (41  ini/hr)  can 
extend  the  range  by  about  1800  ft. 

The  fact  that  a  separate  and  complete  trajectory  is  calculated  for  each  fragment 
permits  the  FRAGHAZ  Computer  Program  to  be  used  for  a  large  variety  of  specific 
problems  like  the  barricade  studies  contained  herein.  The  program  has  also  been 
used  to  determine  hazards  to  vehicles  on  public  traffic  routes  and  hazards  to 
ammunition  stacks  from  fragments. 


MONTE  CARLO  PROCEDURE 


A  Monte  Carlo  procedure  is  any  mathematical  technique  that  uses  random  num¬ 
bers  to  simulate  the  uncertainty  associated  with  one  or  more  random  variables  affect 
ingthe  outcome  of  the  calculations.  In  FRAGHAZ,  there  are  seven  random  variables 
affecting  the  hazard  calculations: 

•  Initial  fragment  elevation  angle 

•  Initial  fragment  velocity 

•  Fragment  drag  coefficient  curve3 

•  Height  of  fragment  trajectory  origin  above  the  ground  surface 

•  Soil  constant  for  ricochet 

•  Windspeed 

•  Altitude  of  the  ammunition  stack  site 

In  the  Monte  Carlo  procedure,  each  replication  represents  one  simulation  of  a  full- 
scale  test.  For  example,  the  trajectories  of  the  260  fragments  used  to  represent  Mk  82 
bomb  pallets  would  constitute  one  replication.  Each  random  variable  associated  with 
the  fragments  would  have  a  known  or  assumed  range  of  uncertainty.  Random  num¬ 
bers  are  then  used  to  designate  a  particular  value  for  each  random  variable.  Any  of 
the  random  variables  may  be  set  to  a  constant  in  FRAGHAZ.  Trajectories  would  be 
calculated  for  each  of  the  260  fragments  with  an  effective  number  of  fragments  asso¬ 
ciated  with  each  trajectory  commensurate  with  the  number  of  munitions  or  inter¬ 
action  areas  in  the  specific  stack  under  consideration.  Hazardous  intersections  with 
the  target  would  be  recorded  and  accumulated  in  the  program  for  each  100-ft  sector  of 
the  hazard  volume.  This  would  constitute  one  replication.  Because  of  the 
uncertainty  in  the  random  variables,  this  would  constitute  only  one  possible  outcome 
for  the  stack  being  considered.  As  a  result,  a  second  replication  would  be  conducted 
with  the  260  fragments,  using  a  new  set  of  random  numbers  to  define  new  values  for 
the  random  variables.  New  outcomes  for  the  100-ft  sectors  of  the  hazard  volume 
would  be  obtained  and  recorded.  This  procedure  would  continue  for  the  number  of 
replications  selected,  60  for  the  data  in  this  report.  When  all  the  replications  are 
completed,  we  are  then  in  a  position  to  calculate  statistics,  such  as  averages  for  such 
measures  of  effectiveness  as  hazard  density  for  each  100-ft  range  increment  of  the 
hazard  volume. 
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For  this  report,  percentiles  are  important.  In  our  example,  we  would  have  60 
distinct  hazard  densities  for  each  100-ft  range  increment.  For  each  100-ft  range 
increment,  we  can  sort  the  60  values  in  ascending  order  from  minimum  to  maximum. 
When  sorted,  the  54th  highest  value  would  be  the  90th  percentile  value  because  only 
6  values  (10  percent  of  the  total)  would  be  equal  to  or  greater  than  the  54th  value. 
The  100th  percentile  would  then  coincide  with  the  largest  sorted  value.  We  may 
interpret  the  percentile  level  as  follows:  If  we  are  talking  about  a  90th  percentile  val¬ 
ue,  we  may  say  that  only  10  percent  of  the  time  should  we  expect  the  value  to  be 
greater  than  the  90th  percentile  value,  and  we  should  not  expect  the  value  to  be 
greater  than  the  100th  percentile  value. 


HAZARD  CRITERIA 


In  order  to  distinguish  hazardous  from  nonhazardous  fragments,  hazard  criteria 
must  be  specified  beforehand.  The  following  criteria  represent  the  current 
specifications  of  DDESB: 


Kinetic  Energy 


A  hazardous  fragment  is  one  that  has  at  least  58  ft-lb  of  kinetic  energy  when  it 
strikes  a  personnel  target.  DDESB  is  currently  reviewing  a  criterion  submitted  by 
the  Army  that  depends  not  only  on  mass  and  velocity  but  on  the  average  presented 
area  of  the  fragment  on  impact.  This  criterion.  Continuous  Probability  of  Injury  Cri¬ 
terion  (CPIC),  which  depends  on  skin  penetration,  has  a  great  deal  of  experiment,.*! 
data  to  support  it.  Some  recent  studies,  however,  indicate  that  the  two  criteria 
produce  roughly  the  same  results. 


Hazard  Density 


Densities  equal  to  or  greater  than  one  fragment  per  600  ft2  are  considered  hazard¬ 
ous,  provided  all  the  fragments  involved  are  hazardous  in  terms  of  the  kinetic  energy 
criterion  given  above.  This  density  is  predicated  on  a  standing  man  having  a  frontal 
presented  area  of  6  ft2.  As  such,  it  is  roughly  equivalent  to  a  probability  of  hazard  hit 
of  0.01. 
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Other 


Both  wind  and  percentile  level  affect  hazard  density  statistics.  Currently, 
DDESB  specifies  zero  wind  and  a  90th  percentile  levels  for  safety  studies. 


PROGRAM  OUTPUT 


The  FRAG11AZ  Computer  Program  has  three  basic  outputs:  Number  of  final 
ground  impacts  versus  range,  hazard  density  and  probability  of  hit  versus  range,  and 
number  of  units  required  to  exceed  the  density  and  probability  of  hit  hazard  criteria  1 
versus  range. 


Number  of  Final  Ground  Impacts  Versus  Range 


If  we  have  60  replications  of  260  fragments,  then  we  will  obtain  60  distinct  distri¬ 
butions  of  final  ground  impacts  versus  range.  This  data  can  be  sorted;  and,  from  it, 
we  can  obtain  minimum,  maximum,  and  average  numbers  of  final  ground  impacts  for 
each  100-ft  increment  of  range.  This  data  can  be  used  to  demonstrate  the  validity  of 
the  program,  as  explained  below  under  the  section,  FRAGHAZ  COMPUTER 
PROGRAM  VALIDATION. 


Hazard  Density  and  Probability  of  Hit 


Hazard  density  is  obtained  by  projecting  the  particular  three  dimensional,  100-ft 
hazard  sector,  intersected  by  the  fragment  trajectory,  into  the  plane  perpendicular  to 
the  trajectory  at  the  point  of  target  impact.  This  can  be  done  since  we  know  the  tra¬ 
jectory  angle  at  all  times.  The  projected  target  area  is  not  calculated  but  is  assumed 
to  be  a  constant  6  ft2  for  a  personnel  target.  The  hazard  density  is  then  the  number  of 
hazardous  fragments  associated  with  the  particular  trajectory  divided  by  the  project¬ 
ed  area  of  the  100-ft  hazard  sector  in  the  plane  perpendicular  to  the  fragment 
trajectory. 

For  probability  of  a  hazardous  hit,  we  use  the  hazard  density  calculated  above  and 
the  presented  area  of  the  target  in  the  plane  perpendicular  to  the  fragment  trajectory. 
These  two  values  are  then  used,  with  an  appropriate  equation,  to  calculate  hazard 
probability  of  hit  for  the  personnel  target. 

Because  target-presented  area  is  not  used  in  the  calculation  of  hazard  density,  us¬ 
ing  the  hazard  density  criterion  will  make  the  results  look  more  hazardous  than  those 
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calculated  using  a  consistent  probability  of  hazard  hit  criterion.  This  occurs  because 
the  hazard  density  criterion  of  one  fragment  per  600  ft'  assumes  the  maximum  pre¬ 
sented  area  of  a  personnel  target,  6  ft’,  which  approximates  a  probability  of  hit  of 
0.01 .  With  a  probability  of  hit  criterion,  the  presented  area  of  the  personnel  target  is 
taken  into  account  and  will  almost  always  be  less  than  6  ft2.  In  fact,  it  can  be  as  low 
as  0.55  ft2  when  the  fragment  comes  straight  down  on  a  standing  man.  As  such,  with 
a  hazard  density  of  1/600,  the  probability  of  hit  will  almost  always  be  less  than  0.01. 


Number  of  Units  To  Exceed  the  Hazard  Criterion 


Using  the  hazard  density  and  probability  of  hit  discussed  above,  one  can  use  the 
program  to  calculate  the  number  of  units  or  interaction  areas  required  to  just  exceed 
the  hazard  density  criterion  (1/600)  or  the  hazard  probability  of  hit  criterion  (0.01). 


FRAGHAZ  COMPUTER  PROGRAM  VALIDATION 


With  computer  programs,  there  is  always  a  question  of  whether  the  predicted 
results  obtained  with  the  program  actually  match  the  results  obtained  with  experi¬ 
mental  tests.  Two  tests  have  been  conducted,  one  with  36  pallets  (288  rounds)  of 
155mm  projectiles  and  another  with  1  pallet  (6  bombs)  of  Mk  82  bombs.  After  the 
explosion  in  each  case,  personnel  were  sent  down  range  to  pick  up  all  fragments 
greater  than  300  grains.  The  terrain  was  marked  off  in  10-deg  azimuthal  sectors;  and 
these  sectors,  in  turn,  were  marked  off  into  200-  to  400-ft  range  increments.  In  this 
way,  each  fragment  could  be  identified  with  a  particular  azimuth  sector  and  a  range 
zone  in  the  azimuth  sector. 

As  mentioned  above,  the  FRAGH  AZ  Computer  Program  can  output  the  number  of 
final  fragment  impacts  as  a  function  of  100-ft  range  increments  for  a  particular  azi¬ 
muth  sector.  Since  we  had  only  one  test  for  155mm  projectiles  and  six  for  Mk  82 
bombs,  we  cannot  say  whether  the  results  per  range  increment  were  maximum  or 
minimum  values  or  some  point  in  between.  Since  each  replication  of  FRAGHAZ  rep¬ 
resents  the  results  of  an  entire  stack-and  we  run  many  replications  (usually  60)--, we 
can  obtain  the  minimum  and  maximum  results  expected  in  the  actual  tests.  We, 
therefore,  ask  that  the  actual  test  results  fall  between  the  minimum  and  maximum 
values  predicted  by  the  FRAGHAZ  Computer  Program.  In  both  cases  there  was  good 
agreement  between  the  predicted  and  actual  test  results. 
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SUMMARY 


The  FRAGHAZ  Computer  Program  provides  a  flexible  tool  for  predicting  the  frag¬ 
ment  hazards  of  stacks  of  ammunition.  The  program  has  the  inherent  capability  of 
considering  the  multidimensional  problem  posed  by  fragment  hazards.  The  program 
is  moderately  large,  having  about  1000  lines  of  code  and  over  200  variables,  about  25 
of  which  are  prime  variables  directly  affecting  hazard  values.  Its  modular  character¬ 
istics  make  it  relatively  easy  to  modify  for  specific  problems  like  the  barricades 
considered  herein.  The  essential  characteristics  of  the  program  are  summarized  as 
follows: 

•  Individual  three-dimensional  trajectories 

•  Two-dimensional  wind  vectors  (horizontal  plane) 

•  Fourth-Order  Runge-Kutta  trajectory  calculations 

•  Fragment  ricochet  included  for  various  soil  types 

•  Incorporates  three-dimensional  targets 

•  Can  use  different  hazard  criteria 

•  Air  density  and  sound  speed,  a  function  of  altitude 

•  Storage  sites  that  may  be  at  different  altitudes 

•  Fragment  drag  coefficient,  a  function  of  Mach  Number  and  based  on  wind 
tunnel  tests 

•  Predicts  distribution  of  final  fragment  impacts  in  ground  plane 

•  Predicts  hazard  densitj'  and  hazard  probability  of  hit  as  a  function  of  range  for 
different  hazard  levels  such  as  minimum,  maximum,  average,  and  specified 
percentiles 

•  Predicts  hazard  distance  values  for  different  hazard  levels  (minimum,  maxi¬ 
mum,  average,  and  specified  percentiles)  as  a  function  of  number  of  units  or 
interaction  areas  required  to  just  exceed  each  of  the  two  hazard  criteria, 
density,  and  probability  of  hit 
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PART  I-GENERAL  BARRICADE  INVESTIGATION 


PREFACE 


This  study  was  requested  by  DDESB,  who  wanted  to  know  if  the  FRAGHAZ 
Computer  Program  could  be  used  to  assist  in  the  design  and  evaluation  of  barricades 
intended  for  use  in  stopping  fragmentation  produced  by  ammunition  stacks.  In  the 
past,  barricade  studies  were  mostly  qualitative  and  based  on  the  extensive 
experience  of  investigators  in  general  fragmentation  characteristics. 

A  related  purpose  of  the  study  was  to  demonstrate  that  meaningful  measures  of 
effectiveness  could  be  produced  by  the  FRAGHAZ  Program  to  assist  in  the  design  and 
evaluation  of  barricades.  The  FRAGHAZ  Computer  Program  was  modified  to  include 
the  geometry  and  variables  associated  with  the  barricade  defined  below. 


CONDITIONS 


Barricade  Model 


The  barricade  model  is  shown  in  Figure  1.  In  keeping  with  the  general  nature  of 
the  study,  the  barricade  used  was  a  simple  box  type.  The  walls  of  the  box  barricade 
could  be  planks,  the  inside  of  which  could  be  filled  with  dirt.  The  variables  addressed 
in  the  study  are  also  shown  in  Figure  1.  The  Stack  Inert  Ground  Standoff  (SIGS)  is 
just  the  height  of  the  pallet  on  which  the  munitions  rest.  The  Barricade  Face  Angle 
(BFA)  determines  the  slope  of  the  front  face  of  the  barricade.  For  this  general  study, 
the  front  face  of  the  barricade  was  vertical  and  the  barricade  face  angle  was  90  deg. 
All  fragments  that  strike  the  barricade  face  or  the  ground  between  stack  and 
barricade  are  considered  stopped  in  the  0-  to  100-ft  range,  increment  of  the  FRAGHAZ 
hazard  volume.  The  barricade  face  angle  is  too  large  to  permit  ricochet;  therefore, 
ricochet  off  the  barricade  was  ignored.  Under  instruction  of  DDESB,  the  complete  or 
partial  destruction  of  the  barricade  by  the  explosion  of  the  munitions  was  not 
considered. 
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HS  -  Height  ot  Stack 
SIGS  -  Stack  Inert  Ground  Standotl 
BS  -  Barricade  Standoff 
BIA  -  Barricade  Intercept  Angle 
BFA  -  Barricade  Face  Angle  (90  deg) 
HB  -  Height  of  Barricade 


FIGURE  1.  BOX  BARRICADE 


Munitions  Selection 


The  two  munitions  selected  were  those  for  which  there  are  fragmentation  data 
appropriate  to  the  FRAGHAZ  Computer  Program.  These  munitions  were  the  Mk  82 
bomb  and  the  M107  155mm  projectile.  Of  all  the  fragmentation  data  available,  these 
two  munitions  are  the  most  hazardous.  Additionally,  the  Mk  82  bombs  are  stored 
horizontally,  while  the  155mm  projectiles  are  stored  vertically.  The  fragmentation 
data  available  reflect  this  difference  in  storage  attitude.  Both  sets  of  fragmentation 
data  are  for  mass  detonating  munitions  (Class  1,  Division  1).  For  both  munitions, 
only  one  pallet  was  used.  For  both  munitions,  the  stack  height  is  approximately  2.5 
ft,  including  the  pallet.  As  such,  the  Barricade  Height  is  determined  by  this  height 
and  the  values  of  Barricade  Intercept  Angle  (BIA)  and  Barricade  Standoff  (BS),  as 
shown  in  Figure  1. 


Hazard  Criteria 


Throughout  this  study,  the  target  was  a  standing  man.  The  CPIC  criterion  was 
used.  This  criterion  depends  on  skin  penetration  and  was  provided  by  the  Ballistic 
Research  Laboratories  at  Aberdeen,  Maryland.  It  is  currently  under  review  by 
DDESB,  Unlike  the  58  ft-lb  criterion,  with  the  CPIC  criterion,  all  fragments  are 
hazardous  to  some  finite  level  of  probability;  therefore,  the  maximum  range  and  the 
maximum  hazard  range  will  be  the  same.  This  is  an  advantage  for  the  measures  of 
effectiveness  selected,  as  described  below. 

The  current  DDESB  hazard  density  criterion  of  at  least  one  hazardous  fragment 
per  600  ft2  is  used.  In  addition,  zero  wind  is  used  along  with  90th  percentile  levels. 
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Measures  of  Effectiveness 


The  two  measures  of  effectiveness  and  associated  criteria  used  are  described 
below: 

•  The  first  was  maximum  hazard  range.  Since  we  are  using  the  CPIC  injury 
Criterion,  all  fragments  are  hazardous,  and  the  maximum  fragment  range 
coincides  with  the  maximum  hazard  range.  This  is  unlike  the  case  with  the  58 
ft-lb  injury  criterion,  where  the  fragment  going  to  maximum  range  may  have  a 
kinetic  energy  less  than  58  ft-lb  and  therefore  would  not  be  considered  hazard¬ 
ous  by  the  program.  As  such,  maximum  range  and  maximum  hazard  range 
could  be  different. 

•  The  second  measure  of  effectiveness  is  hazard  density.  This  is  the  measure 
currently  specified  by  DD/SSB.  For  a  100-ft  range  increment  to  be  hazardous, 
the  density  must  be  at  least  one  fragment  per  600  ft2  (1/600  =  0.001667),  and 
all  fragments  making  up  the  density  must  be  hazardous.  Since  we  are  using 
the  CPIC  injury  criterion,  all  fragments  are  hazardous  to  a  finite  probability, 
and  this  probability  is  used  in  calculating  the  reduced,  or  effective,  hazard 
density. 

•  The  two  measures  of  effectiveness  were  calculated  using  zero  wind  and  90th 
percentile  levels.  The  barricade  standoff  (Figure  1)  was  held  constant  at  4  ft, 
the  practical  minimum  specified  by  DDESB.  Hazards  will  tend  to  increase  as 
BS  increases. 


Monte  Carlo  Running  Option 


The  Monte  Carlo  option  was  used  in  all  cases.  All  runs  contained  60  replications. 
The  same  seed  (starting  number  used  by  the  random  number  generator)  (one  for  each 
munition)  was  used  such  that  the  same  fragment  trajectories  were  used  for  each 
barricade  intercept  angle.  In  this  way,  the  dispersion  in  results  with  different  seeds 
was  eliminated;  and  only  the  effects  of  different  barricade  intercept  angles  are  seen. 


Stack  Top  Fragmentation 


Because  of  the  physical  dimensions  of  the  munition  pallets,  a  part  of  the  frag¬ 
mentation  from  the  tops  of  the  stacks  was  included  in  the  arena  test  recovery.  The 
exact  amount  of  the  fragmentation  from  the  stack  tops,  however,  is  unknown.  For 
this  initial  study,  it  is  not  of  much  consequence,  since  the  elevation  angles  involved 
would  be  greater  than  70  deg.  In  Part  n,  however,  where  stacks  may  be  6  pallets 
deep,  stack  top  fragmentation  must  be  taken  into  account. 
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RESULTS 


Maximum  Hazard  Range 


The  maximum  hazard  range  as  a  function  of  Barricade  Intercept  Angle  is  shown 
in  Figures  2  and  3  for  155mm  projectiles  and  Mk  82  bombs,  respectively.  The  data  on 
these  two  figures  represent  the  90th  percentile  level.  Maximum  hazard  ranges  for 
the  no-barricade  condition  are  also  shown  for  reference.  Barricade  heights  are  also 
shown  for  corresponding  barricade  intercept  angles,  given  a  stack  height  of  2.5  ft.  If 
the  stack  height  were  increased,  say  to  3  pallets  high  or  7.5  ft,  then  the  barricade 
heights  would  increase  5  ft  for  the  given  barricade  intercept  angles. 


Barricade  Intercept  Angle  (BIA)  •  deg 


| — ' | — ■ — I — ■ 1 — r~r-i — I — ■ I ’ 1 — •— 1 

2.5  3.2  4.0  4.0  5.8  7.3  9.4  13.5  25.2 

Barricade  Height  (HB)  - It 

FIGURE  2.  BARRICADE  EFFECTIVENESS  M107  PROJECTILES  (155mm) 
MASS  DETONATION  (CLASS  1.  DIVISION  1)  VERTICAL  STORAGE 


The  figures  show  a  slow  and  steady  decline  in  maximum  hazard  range  for  in¬ 
creasing  Barricade  Intercept  Angle.  There  are  no  abrupt  changes  except  at  0-deg 
Barricade  Intercept  Angle,  where  we  experience  drops  from  the  no-barricade 
conditions.  At  these  conditions,  the  barricade  height  is  2.5  ft  and  blocks  the  long- 
range  fragments.  In  actuality,  the  two  curves  will  meet  to  the  left  of  the  ordinate  at  a 
barricade  intercept  angle  of  about  —32  deg.  Negative  angles  were  not  considered 
because  of  the  need  to  take  ricochet  into  account  from  the  top  of  the  barricade. 

Note  in  Figures  2  and  3  that  a  small  Barricade  Intercept  Angle  (like  the  2  deg  in 
DOD  6055.9-STD)  does  not  produce  a  drastic  reduction  in  maximum  hazard  range  as 
might  have  been  expected.  The  drops,  however,  are  significant,  about  500  ft  for  the 
155mm  projectiles  and  950  ft  for  the  Mk  82  bombs. 
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Bwricad*  Iniaroapt  Angto  (RtA)  -  (teg 
i — ■  i — » — i — • — > — » — i — • — i — 1 — r — • — i — • — i 
2.5  3.2  4.0  4.8  5.9  7.3  9.4  13.5  25.2 

Barricade  Height  (HB)  -  it 

FIGURE  3.  BARRICADE  EFFECTIVENESS  MK  82  LOW  DRAG  BOMBS  (500  lb) 
MASS  DETONATION  (CLASS  1,  DIVISION  1)  HORIZONTAL  STORAGE 


Table  1  gives  maximum  hazard  ranges  for  the  50th,  90th,  and  100th  percentile 
levels. 

Maximum  hazard  range  provides  a  simple  and  efficient  measure  of  barricade 
effectiveness.  It  does  not  depend  on  the  number  of  munitions  in  the  stack  as  hazard 
density  does. 


Hazard  Density 


Hazard  density,  as  a  function  of  range  and  Barricade  Intercept  Angle,  is  shown  in 
Figures  4  and  5  for  155mm  projectiles  and  Mk  82  bombs,  respectively.  Again,  the 
data  in  these  figures  represent  90th  percentile  levels.  The  current  DDESB  density 
criterion  (1  frag/600  ft2)  is  given  on  the  figures  for  reference.  Hazard  density  depeiids 
on  the  number  of  munitions  or  interaction  areas  on  the  face  of  the  stack  toward  the 
target  area.  The  data  shown  on  Figures  4  and  5  are  for  one  interaction  area.  If  haz¬ 
ard  densities  for  more  than  one  interaction  area  are  desired,  the  values  given  in 
Figures  4  and  5  can  be  multiplied  by  the  number  of  interaction  areas  in  question.  For 
example,  in  Figure  5  the  hazard  density  at  a  range  of  2000  ft  and  a  BIA  of  0  deg  is 
approximately  0.00001.  If  a  stack  had  50  interaction  areas  on  the  face  of  the  stack  to¬ 
ward  the  target  area,  then  the  hazard  density  would  be  0.0005  at  a  range  of  2000  ft 
and  a  BIA  of  0  deg.  If  we  wished  to  know  how  many  interaction  areas  on  the  face  of 
the  stack  toward  the  target  area  would  result  in  a  hazard  density  just  greater  than 
the  DDESB  criterion  of  1  frag/600  ft2  (0.001667),  we  would  divide  the  hazard  criterion 
by  the  hazard  density  for  one  interaction  area.  For  the  case  stated  of  a  BIA  of  0  deg  at 
a  range  of  2000  ft,  this  would  be  (0.001667/0.00001)  or  167  when  rounded  up  to  the 


13 


NAVSWC  MP  89-353 


TABLE  1.  MAXIMUM  HA7-AHI)  RANGE  (ft)  FOR  BOX  BARRICADES* 


Barricade 
Intercept  Angle 
(deg) 

Munition:  155  mm  Projectile 

Percentile: 

50  90  100 

Munition:  Mk  82  Bomb 
Percentile: 

50  90  100 

Barricade 
Height  (ft) 

0 

1900 

2400 

2400 

1900 

2400 

2400 

2.5 

5 

1500 

2400 

1900 

2400 

2400 

29 

10 

1400 

Bs 

2400 

1800 

2300 

2400 

3.2 

20 

1300 

1600 

2400 

1700 

2000 

2400 

4.0 

30 

1200 

1400 

1600 

1500 

1800 

2000 

4.8 

40 

1200 

1400 

1600 

1200 

1700 

1800 

5.9 

50 

1100 

1400 

1600 

1100 

1300 

1700 

7.3 

60 

900 

1200 

1500 

900 

1200 

1400 

9.4 

70 

700 

900 

1000 

700 

800 

1000 

13.5 

80 

400 

600 

600 

300 

300 

400 

25.2 

No  barricade 

2300 

2900 

3200 

2000 

3300 

3500 

•Conditions  for  both  munitions: 

Height  of  stack  =  2.5ft 
Stack  inert  ground  standoff  =  0.6  ft 
Barricade  standoff  =  4.0  ft 
Barricade  face  angle  =  90  deg 
Barricade  intercept  angle  =  0-80  deg 

nearest  whole  number.  Without  the  barricade,  this  number  would  be  reduced  to 
about  (0.001667/0.00003)  or  56  interaction  areas,  which  we  might  say  is  about  three 
times  more  hazardous. 

Hazard  density  might  be  more  properly  called  injury  density  when  using  the 
CPIC.  The  hazard  density  criterion  (1/600)  is  used  to  approximate  the  probability  of 
hit  criterion.  The  hazard  density  is  defined  as  follows: 

Hazard  density  =  (NF/AT)  *  P  (INJ/HIT) 

where 


NF  =  Number  of  hazardous  fragments 

AT  =  Presented  area  of  the  hazard  volume  of  a  100-ft  range  segment  in 
the  plane  perpendicular  to  the  fragment  trajectory 
P  (INJ/HIT)  =  Conditional  probability  of  injury  given  a  hit 
*  =  Multiplication 
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FIGURE  4.  BARRICADE  EFFECTIVENESS  M107  PROJECTILES  (155mm) 
MASS  DETONATION  (CLASS  1,  DIVISION  1)  VERTICAL  STORAGE 


Hazard  Rang*  (R) 


FIGURE  5.  BARRICADE  EFFECTIVENESS  MK  82  LOW-DRAG  BOMBS  (500  lb) 
MASS  DETONATION  (CLASS  1,  DIVISION  1)  HORIZONTAL  STORAGE 
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With  the  58  ft-lb  injury  criterion,  the  probability  of  injury  is  either  zero  or  one  as 
the  impacting  fragment  has  less  than  58  ft-lb  or  at  least  58  ft-lb  of  kinetic  energy, 
respectively. 


CONCLUSIONS 


The  FRAGHAZ  Computer  Program  can  be  adapted  to  evaluate  the  design  and 
effectiveness  of  barricades. 

More  complex  barricade  designs  and  effects  of  ricochet  can  be  considered  by  the 
FRAGHAZ  Computer  Program. 

The  FRAGHAZ  Computer  Program  can  be  used  to  determine  shapes  and  dimen¬ 
sions  for  barricades  to  meet  specific  design  criteria.  For  example,  it  might  be  desir¬ 
able  to  design  a  barricade  to  limit  down-range  hazardous  fragments  to  the  hazard 
blast  radius  for  a  specific  ammunition  stack. 
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PART  IE-SPECIFIC  HILLSIDE  BARRICADE  INVESTIGATION 


PREFACE 


The  investigation  contained  herein  was  conducted  in  response  to  a  request  made 
by  the  Air  Force.  The  letter  containing  the  request  is  reproduced  in  the  appendix. 
Correspondence  subsequent  to  the  request  provided  details  on  the  site  necessary  to 
conduct  the  investigation.  The  site  is  located  at  Taegu  Air  Base  in  Korea. 

The  Air  Force  requested  that  the  clear  zone  for  fragments  about  a  storage  shed  be 
reduced  from  1250  ft  to  something  nearer  800  ft  because  of  a  natural  hill  barricade 
surrounding  the  shed  on  three  sides.  Although  the  shed  is  used  for  handling  a  variety 
of  munitions,  it  has  been  stated  that  six  Mk  82  bomb  pallets  (36  bombs)  would  be 
taken  as  the  source  of  the  maximum  hazard  for  the  investigation. 

The  FRAGHAZ  Computer  Program  was  modified  to  include  the  site  geometry  and 
munition  fragmentation  characteristics  of  Mk  82  bomb  pallets.  The  Mk  82  bomb 
weighs  500  lb. 


CONDITIONS 


Stack  Configurations 


Six  stack  configurations  were  used,  as  shown  in  Figure  6.  They  were  selected  to 
provide  minimum  and  maximum  fragmentation  hazards  for  both  side  and  top  inter¬ 
action  areas.  Configurations  1  and  5  provide  maximum  and  minimum  hazards  for 
side  interaction  areas,  respectively.  Configurations  3  and  5  provide  minimum  and 
maximum  hazards  for  top  interaction  areas,  respectively.  The  remaining  config¬ 
urations  provide  intermediate  hazards  that  give  information  on  transition  from 
minimum  to  maximum  hazards. 

Information  from  the  Air  Force  stated  that  a  maximum  of  six  Mk  82  bomb  pallets 
(36  bombs)  could  be  stored  in  the  shed  (Bldg.  230)  at  Taegu  Air  Base  in  Korea.  The 
six  configurations,  shown  in  Figure  6,  represent  rectangular  parallelepiped  stacking 
of  the  6  pallets.  Configurations  where  the  number  of  pallets  would  differ  in  different 
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FIGURE  6.  STACK  CONFIGURATIONS,  KOREAN  HILLSIDE  BARRICADE,  TAEGU,  BLDG.  230, 

6  PALLETS  Mk  82,  6  BOMBS  PER  PALLET 


NAVSWC  MP  89-353 


tiers  were  not  considered,  since  they  would  only  involve  intermediate  hazard 
conditions  between  the  six  configurations  selected. 

Because  the  height  of  the  shed  is  about  10  ft,  the  maximum  height  of  any  stack 
was  taken  to  be  7.5  ft;  that  is,  3  pallets  high.  The  height  of  an  individual  pallet  is 
about  2.5  ft. 

The  abbreviations  associated  with  each  stack  configuration  in  Figure  6  are 
explained  as  follows: 

HS  =  Height  of  stack 

SIA  —  Number  of  side  interaction  areas  used  in  the  calculation  of  hazard  den¬ 
sities  (With  the  hill  barricade  (w/hill),  only  the  interaction  areas  of  the 
top  layer  of  the  pallets  are  used.  The  fragmentation  from  the  inter¬ 
action  areas  of  lower  layers  will  be  stopped  by  the  hillside  barricade. 
The  top  layer  includes  not  only  the  interaction  area  fragments  from  the 
bomb  case  zone  A-B-C  in  Figure  7  but  also  fragmentation  from  the  zone 
C-D.  Without  the  hill  (w/o  hill),  all  interaction  areas  on  the  face  of  the 
stack,  including  those  formed  by  adjacent  pallets,  are  pertinent  to  the 
calculations.  When  pallets  are  stacked  on  top  of  one  another,  an 
additional  interaction  area  is  formed  between  pallets.) 

TLA  ~  Number  of  interaction  areas  on  the  top  of  the  stack,  including  those 
formed  by  adjacent  pallets  (When  pallets  are  positioned  next  to  one 
another  with  long  bomb  axis  to  long  bomb  axis,  then  an  additional  inter¬ 
action  area  is  formed  between  pallets.  If  the  pallets  are  positioned  bomb 
base  to  bomb  base,  nose  to  nose,  or  nose  to  base,  then  no  additional 
interaction  areas  are  formed.) 

The  effects  of  blast  pressures  from  any  of  the  stacks  were  not  considered  in  the 
hazard  calculations. 


Site  Configuration 


The  essentials  of  the  site  configuration  used  for  hazard  calculations  are  shown  in 
Figure  6.  The  hillside  barricade  consists  of  two  legs,  one  vertical  and  one  angled  back 
at  60  deg.  Fragment  ricochet  was  considered  on  the  leg  angled  back  at  60  deg.  A  soil 
constant  of  l.l  was  used,  which  is  appropriate  for  cohesive  soils  supporting  vegetative 
growth. 

The  effects  of  the  shed  structure  on  fragment  trajectories  were  not  considered. 
The  effects  are  considered  minor.  Ignoring  shed  effects  might  even  tend  to  err 
against  the  side  of  safety,  since  the  shed  structure  could  be  the  source  of  additional 
fragments,  although  such  fragments  would  probably  have  poor  aerodynamic 
characteristics. 
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FIGURE  7.  FRAGMENTATION  ELEVATION  ZONES  (KOREAN 
HILLSIDE  BARRICADE,  TAEGU,  BLDG.  230) 


The  site  altitude  used  in  the  hazard  calculations  was  150  ft  above  sea  level-  Even 
if  this  altitude  is  off  by  a  few  hundred  feet,  the  effects  on  the  overall  calculations 
would  be  negligible. 

Three  distances  for  positioning  the  stacks  relative  to  the  hillside  were  used.  A  BS 
of  32  ft,  as  shown  in  Figure  6,  corresponds  to  a  stack  position  at  the  back  of  the  shed. 
A  BS  of  15  ft  is  for  a  stack  position  at  the  front  of  the  shed.  A  BS  of  5  ft  was  used  to 
determine  whether  a  significant  reduction  in  fragment  hazards  could  be  obtained 
that  might  justify  a  change  in  storage  location. 


Hazard  Criteria 


The  target  for  this  investigation  was  a  standing  man.  The  58  ft-lb  injury  cri¬ 
terion,  currently  specified  by  DDESB,  was  used  to  distinguish  hazardous  from  non- 
hazardous  fragments.  The  CPIC  criterion  was  used  for  comparison  purposes.  The 
CPIC  criterion  was  provided  by  the  Army  and  is  currently  under  review  by  DDESB. 

A  hazard  density  of  one  hazardous  fragment  per  600  ft2,  currently  specified  by 
DDESB,  was  used  as  the  limit  above  that  which  hazard  density  is  considered  unac¬ 
ceptable  for  safety  purposes  by  DDESB. 

Zero  wind  and  the  90th  percentile  level  were  also  used  as  standard  conditions  for 
safety  analysis  as  currently  specified  by  DDESB. 
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Side  and  Top  Fragmentation 


The  zones  for  side  and  top  fragmentation  are  shown  in  Figure  7.  Side  frag¬ 
mentation  comes  from  the  bomb  case  portion  labeled  Z-A-B-C-D.  The  side  interaction 
area  is  taken  to  come  from  A-B-C,  where  the  fragment  densities  and  velocities  ore 
highest.  In  the  arena  test,  from  which  the  fragmentation  characteristics  used  in  this 
investigation  were  tak^n,  the  pallet  of  Mk  82  bombs  was  positioned  so  as  to  recover 
the  fragmentation  from  the  side  area.  At  the  time,  it  was  assumed  that  all  fragmen¬ 
tation  going  down  range  came  from  the  face  of  the  stack  toward  the  target  area. 
Reviewing  the  arena  test  data,  however,  showed  that  fragments  with  interaction 
area  velocities  were  recovered  between  polar  angles  10  and  30  deg.  These  fragments 
must  almost  certainly  come  from  the  top  interaction  area  zone  D-E-F.  Whether  some 
of  this  fragmentation  might  also  have  come  from  top  interaction  area  zone  F-G-H  is 
unknown.  Because  of  this  uncertainty  and  for  the  sake  of  safety,  the  side  and  top 
fragmentations  were  taken  to  be  exclusive  of  one  another.  That  is,  the  fragmentation 
recovered  in  the  arena  test  was  taken  to  come  from  the  bomb  case  zone  Z-A-B-C-I) 
only.  Adopting  such  a  procedure  may  result  in  some  double  counting  of  fragments, 
but  any  error  should  tend  to  be  on  the  side  of  safety. 

The  top  fragment  interaction  areas  were  taken  to  be  similar  to  the  side  interaction 
area  A-B-C.  In  reality,  much  depends  on  where  the  stack  is  initiated.  In  the  arena 
test,  the  middle  bomb  of  the  lower  tier  was  the  donor  bomb  making  the  interaction 
area  A-B-C  to  be  skewed  upwards.  In  practice,  the  initiation  point  for  a  stack  will  be 
unknown.  As  such,  we  have  used  an  average,  such  that  the  bomb  case  zone  A-B-C  is 
taken  to  produce  an  interaction  area  between  polar  angles  70  and  110  deg  (elevation 
angles  plus  or  minus  20  deg).  This  zone  then  becomes  the  model  for  the  top  interac¬ 
tion  areas.  Each  top  interaction  area  is  taken  to  be  independent,  since  we  have  no 
test  data  on  the  interaction  of  interaction  areas.  The  arena  test  was  conducted  with 
only  one  interaction  area,  A-B-C.  However,  this  same  type  of  independence  was  used 
with  155mm  projectiles,  which  had  six  interaction  areas  in  the  arena  test;  and  good 
results  were  obtained. 

The  top  interaction  areas  were  defined  by  rotating  the  side  interaction  area,  A-B- 
C,  90  deg,  such  that  elevation  angle  -20  deg  became  elevation  angle  70;  and  eleva¬ 
tion  angle  20  deg  became  elevation  angle  110.  The  top  interaction  areas  have  a 
spread  of  40  deg  like  the  side  interaction  area,  When  making  such  a  90-deg  rotation, 
account  must  be  taken  of  a  diminished  number  of  fragments  that  are  applicable  to  the 
10-deg  hazard  volume.  When  recovering  fragments  in  the  arena,  the  vertical  extent 
of  the  recovery  packs  at  any  elevation  angle  between  the  10-deg  azimuthal  limits  is 
approximately 

H  =  2*  R  *  TAN  (5°)  ♦  ABS  (COS  (E)) 
where 

H  =  Distance  between  azimuthal  limits 

R  =  Radius  of  the  arena 

E  =  Elevation  angle  in  degrees 

*  =  Multiplication 
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When  we  rotate  toward  polar  angle  0  (elevation  angle  90),  the  distance  between 
the  azimuthal  limits  decreases,  finally  becoming  zero  at  polar  angle  0  deg  (elevation 
angle  90  deg).  Since  only  those  fragments  whose  initial  trajectories  are  constrained 
between  the  azimuthal  limits  will  end  up  in  the  10-deg  hazard  volume,  we  must  take 
this  into  account  by  diminishing  the  number  of  fragments  recovered  from  the  polar 
zone  70  to  110,  which  we  rotate  90  deg  to  form  top  interaction  areas.  The  reduction 
factors  used  were  0.25  for  the  70-  to  1 10-deg  top  interaction  area  elevation  zone,  and 
0.1  for  the  80  to  90  top  interaction  area  elevation  zone. 

In  the  investigation,  the  hazard  densities  were  calculated  independently  for  the 
side  and  top  fragmentation  and  then  summed  for  total  hazard  density.  Linear  scaling 
was  used  for  multiple  interaction  areas;  that  is,  hazard  densities  were  calculated  for 
one  interaction  area  and  then  multiplied  by  the  number  of  interaction  areas  for  the 
particular  stack  configuration  to  obtain  total  hazard  densities. 


Wind 


DDESB  currently  specifies  zero  wind  for  safety  analysis.  Two  other  wind 
conditions  were  used  for  reference.  These  were  a  30-ft/sec  tailwind  (20.5  mi/hr)  and 
60-ft/sec  tailwind  (40.9  mi/hr).  A  tailwind  tendi  g  to  push  the  fragments  directly 
down  range  is  almost  always  the  most  hazardous.  A  30-ft/sec  wind  corresponds  to  fair 
weather  gusting.  A  60-ft/sec  wind  corresponds  to  gusting  associated  with  electrical 
storms. 


Monte  Carlo  Running  Option 


The  Monte  Carlo  option  of  the  FRAGHAZ  Computer  Program  was  used  in  all 
cases.  All  runs  contained  60  replications.  The  same  seed  was  used  in  all  cases,  such 
that  the  same  fragment  trajectories  were  used  for  all  six  stack  configurations  and  the 
alternate  conditions  studied.  In  this  way,  the  dispersion  in  results  with  different 
seeds  was  eliminated,  and  only  the  effects  of  the  different  stack  configurations  and 
alternate  conditions  were  obtained. 


RESULTS 


In  order  to  provide  a  graphic  picture  of  the  results,  selected  data  have  been  plotted 
in  Figures  8  through  15.  These  figures  show  the  effect  on  results  caused  by  changes 
in  hazard  criteria  and  other  alternate  conditions.  When  plotting  these  figures,  the 
hazard  densities  for  the  0-  to  100-ft  sectors  have  not  been  plotted.  These  densities  are 
extremely  large,  since  they  include  the  hazards  from  all  the  fragments  stopped  by  the 
hill.  A  discussion  of  the  results,  shown  in  Figures  8  through  15,  is  given  below. 
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Baseline  Hazard  Densities  (Figure  8) 


Figure  8  shows  three  curves,  one  for  each  of  the  three  tailwind  speeds  (0,  30,  and 
60  ft/sec).  The  curve  for  a  tailwind  speed  of  0  ft/sec  is  the  baseline  and  will  be  the 
curve  used  as  a  basis  of  comparison  for  most  of  the  succeeding  curves.  The  baseline 
curve  is  considered  the  best  estimate  of  hazard  densities  for  stack  configuration  3  (the 
least  hazardous  configuration)  under  the  hazard  criteria  and  conditions  currently 
specified  by  the  DDESB  as  listed  in  Figure  8.  Note  that  the  barricade  (hill)  standoff" 
is  32  ft,  which  corresponds  to  the  back  of  the  shed.  This  position  causes  a  sightly 
higher  hazard  than  the  other  two  standoffs— 15  and  5  ft.  All  the  remaining  curves  in 
this  figure  and  in  Figures  9  through  15  are  for  comparison  with  this  baseline  curve 
and  to  lend  scope  to  the  investigation  by  showing  the  effects  of  different  hazard 
criteria  and  conditions. 


Target-Standing  Man:  Injury  Criterion-Kinetic  Energy  Threshold:  58ft-lb 
Monte  Carlo  Option:  Replications  60:  Seed:  17351:  Site  Altitude:  150  ft 
Stack  Top  Vertical  Fragmentation  Elevation  Zone:  70  to  110  deg 
90th  Percentile:  Stack  Barricade  (Hill)  Standoff:  32  ft 


0  1000  2000  3000 


1/600 

0*  (Baseline) 
30* 

60* 


*  Tail  Wii.d  Speed  (fl/scc) 


Range  (ft) 

FIGURE  8.  BASELINE  HAZARD  DENSITIES  (KOREAN  HILLSIDE  BARRICADE, 
TAEGU,  BLDG.  230,  STACK  CONFIGURATION  3) 


Note  that  the  ordinate  of  Figure  8  and  succeeding  figures  is  graduated  with  a  log 
scale.  Most  of  us  are  more  accustomed  to  the  usual  linear  scale  graduations,  and  this 
can  be  a  problem  when  we  are  estimating  how  many  times  greater  one  point  is  than 
another.  For  example,  the  threshold  safety  hazard  density  (one  fragment  per  600  ft2), 
shown  as  a  thick  horizontal  line  at  hazard  density  0.001667,  appears  to  be  about  1.3 
times  greater  than  the  baseline  hazard  density  at  200-ft  range.  Actually  the  baseline 
hazard  density  at  200-ft  range  is  0.000449;  and,  therefore,  the  threshold  safety 


23 


NAVSWCMP  89-353 


hazard  density  is  (0.001667/0.000449)  or  3.7  times  greater.  The  threshold  safety 
hazard  density  appears  about  2.  4  times  greater  than  the  baseline  hazard  density  at 
1400-ft  range,  while  actually  it  is  (0.001667/0.000084)  or  19.8  times  greater. 

Since  we  have  always  tended  to  err  on  the  side  of  safety  when  we  encountered  un¬ 
certainties,  the  hazard  densities  for  the  baseline  conditions  represent  a  large  factor  of 
safety  relative  to  the  threshold  safety  hazard  density.  We  might  say  from  the  above 
calculations  that  the  factor  of  safety  ranges  from  about  3  to  20,  depending  upon 
range.  For  the  approximate  safety  range  of  800  ft,  requested  in  the  appendix,  the 
factor  of  safety  would  be  about  18. 

From  Figure  8,  we  can  also  see  that  tailwind  speed  does  not  have  a  very  large  ef¬ 
fect  relative  to  the  threshold  safety  hazard  density.  The  tailwind  curves  for  30  and  60 
ft/sec  are  still  well  below  the  one  fragment  per  600-ft2  (0.001667)  threshold. 


Minimum  and  Maximum  Hazard  Configurations  (Figure  9) 


Figure  9  shows  minimum  and  maximum  hazard  density  curves  for  stack 
configurations  3  and  5,  respectively.  Actually,  there  is  a  crossover  at  about  1200-ft 
range  because  Configuration  5  has  half  as  many  side  interaction  areas  as  Config¬ 
uration  3  (see  discussion  on  Figure  14  for  side  and  top  contributions).  The  Config¬ 
uration  3  curve  is  taken  from  Figure  8.  Configuration  5  (see  Figure  6)  represents  a 
noncompact  stack  of  6  pallets  and  is  somewhat  analogous  to  scattering  single  pallets 
around  the  shed  floor.  Figure  9  indicates  the  need  for  arranging  the  pallets  in 
compact  stacks  like  Configuration  3,  especially  when  hazard  ranges  of  interest  are 
less  than  about  1000  ft. 


Hazard  Densities  for  58  ft-lb  and  CPIC  Injury  Criteria  (Figure  10) 


Figure  10  shows  the  baseline  curve  and  a  second  curve  for  the  same  conditions, 
except  that  the  baseline  personnel  injury  criterion  of  58  ft-lb  has  been  changed  to  the 
CPIC  criterion.  As  shown  in  Figure  10,  the  effect  of  the  change  is  small,  both  curves 
being  well  below  the  threshold  safety  hazard  density  of  one  hazardous  fragment  per 
600  ft2  (0.001667). 


90th  Versus  100th  Percentile  Hazard  Densities  (Figure  11) 


Figure  11  shows  results  with  the  baseline  Configuration  3  stack  at  two  percentile 
levels,  90  and  100.  The  90th  percentile  curve  represents  the  level  of  hazard  densities 
that  we  expect  to  be  exceeded  only  10  percent  of  the  time.  Of  that  10  percent  of  the 
time,  the  100th  percentile  represents  the  maximum  level  expected.  The  curves  are 
close  and  well  below  the  threshold  safety  hazard  density  of  one  hazardous  fragment 
per  600  ft2  (0.001667). 
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Target-Standing  Man:  Injury  Criterion-Kinetic  Energy  Threshold:  58  ft-lb 
Monte  Carlo  Option:  Replications  80:  Seed:  17351:  Site  Altitude:  150  ft 
Stack  Top  Vertical  Fragmentation  Elevation  Zone:  70  to  1 10  deg 
90th  Percentile:  Stack  Barricade  (Hill)  Standoff:  32  ft:  Zero  Wind 


1/600 

Configuration  3  (Baseline) 
Configuration  5 


FIGURE  9.  MINIMUM  AND  MAXIMUM  HAZARD  CONFIGURATIONS 
(KOREAN  HILLSIDE  BARRICADE,  TAEGU,  BLDG.  230) 


Target-Standing  Man 

Monte  Carlo  Option:  Replications  60:  Seed:  17351:  Site  Altitude:  150  ft 
Stack  Top  Vertical  Fragmentation  Elevation  Zone:  70  to  1 10  deg 
90th  Percentile:  Stack  Barricade  (Hill)  Standoff:  32  ft:  Zero  Wind 


Range  (ft) 


1/600 

58  ft  lbs  (Baseline) 
CPIC 


FIGURE  10.  HAZARD  DENSITIES  FOR  58  ft-lb  AND  CPIC  INJURY 
CRITERIA  (KOREAN  HILLSIDE  BARRICADE,  TAEGU,  BLDG.  230, 
STACK  CONFIGURATION  3) 
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Target-Standing  Man:  Injury  Criterion-Kinetic  Energy  Threshold:  58  ft-lb 
Monte  Carlo  Option:  Replications  60:  Seed:  17351:  Site  Altitude:  150  ft 
Stack  Top  Vertical  Fragmentation  Elevation  Zone:  70  to  110  deg 
90th  Percentile:  Stack  Barricade  (Hill)  Standoff:  32  ft:  Zero  Wind 


1/600 

90th  Percentile  (Ba 
100th  Percentile 


Range  (ft) 

FIGURE  11.  90th  VERSUS  100TH  PERCENTILE  HAZARD  DENSITIES  (KOREAN 
HILLSIDE  BARRICADE,  TAEGU,  BLDG.  230,  STACK  CONFIGURATION  3) 


*ne) 


Hazard  Densities  With  and  Without  Hill  (Figure  12) 


Figure  12  shows  the  results  with  stack  Configuration  3  for  the  case  of  no  hill  as  op¬ 
posed  to  the  baseline  case  with  the  hill.  The  case  without  the  hill  implies  a  flat  and 
horizontal  plane  extending  from  the  stack  down  range.  The  curve  for  the  case  of  no 
hill  shows  hazard  density  exceeding  the  DDESB  limit  (one  hazardous  fragment  per 
600  ft2)  out  to  about  1400  ft.  Configuration  3,  however,  is  not  the  best  stack  configu¬ 
ration  for  the  case  without  the  hill.  Another  study  would  have  to  be  performed  to 
define  the  stack  configuration  that  would  be  best  for  open  terrain.  The  use  of 
Configuration  3  is  for  comparison  but  probably  represents  a  near-maximum  hazard 
for  open  terrain. 

In  the  appendix,  the  Air  Force  stated  that  the  storage  shed  is  surrounded  on  three 
sides  by  the  hill.  The  fourth  is  apparently  open  but  blocked  to  some  extent  by  an  un¬ 
defined  barricade.  Since  hazard  densities  for  the  fourth  side  would  probably  fall  be¬ 
tween  the  two  curves  shown  in  Figure  12,  calculations  should  be  made  when  the 
stack-barricade  geometry  for  the  fourth  side  becomes  available. 
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Target-Standing  Man:  Injury  Criterion-Kinetic  Energy  Threshold:  58  ft-lb 
Monte  Carlo  Option:  Replications  60:  Seed:  17351:  Site  Altitude:  150  ft 
Stack  Top  Vertical  Fragmentation  Elevation  Zone:  70  to  110  deg 
90th  Percentile:  Stack  Barricade  (Hill)  Standoff:  33  ft:  Zero  Wind 


1/600 

With  Hill  (Baseline) 
Without  Hill 


FIGURE  12.  HAZARD  DENSITIES  WITH  AND  WITHOUT  HILL  (KOREAN  HILLSIDE 
BARRICADE,  TAEGU,  BLDG.  230,  STACK  CONFIGURATION  3) 


Hazard  Densities  for  0-  and  58-ft-lb  Injury  Criteria  (Figure  13) 


Figure  13  shows  the  effect  of  an  ultra-conservative  personnel  injury  criterion. 
With  the  0-ft-lb  injury  criterion,  every  fragment  has  a  probability  of  injury  of  1. 
Actually,  every  fragment  has  a  kinetic  energy  greater  than  0  at  impact  because  the 
impact  velocity  is  always  greater  than  0.  Selecting  0  ft-lb  as  the  criterion  ensures 
that  every  fragment  will  be  hazardous  when  we  do  not  know  the  minimum  kinetic 
energy  at  impact.  This  figure  is  designed  to  show  an  extreme  point  of  reference. 


Hazard  Density  Contributions  From  Side  and  Top  Fragmentation  (Figure  14) 


It  is  interesting  to  note  that  the  two  contributions  (from  side  and  top)  do  not 
overlay  in  range.  These  two  curves  together  form  the  baseline  curve  but  shows  a 
range  gap.  In  the  previous  figures,  the  baseline  curve  showed  these  two  regions 
connected  to  make  a  smooth,  continuous  curve.  The  figure  shows  that  the  hill  is 
doing  an  excellent  job  of  blocking  the  side  fragmentation.  See  Figure  7  for  a 
description  of  side  and  top  fragmentation. 

If  the  baseline  curve  had  been  calculated  with  the  stack  positioned  at  the  front  of 
the  shed,  the  contribution  from  the  side  fragmentation  would  be  zero  beyond  100-ft 
range.  As  such,  the  baseline  curve  would  terminate  at  900-ft  range. 
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Target-Standing  Man: 

Monte  Carlo  Option:  Replications  60:  Seed:  17351:  Site  Altitude:  150  ft 
Stack  Top  Vertical  Fragmentation  Elevation  Zone:  70  to  110  deg 
90th  Percentile:  Stack  Barricade  (Hill)  Standoff:  32  ft:  Zero  Wind 


1*00 

58  fi-lbs  (Baseline) 
Ofi-lbs 


FIGURE  13.  HAZARD  DENSITIES  FOR  9-  AND  58-ft-lb  INJURY  CRITERIA  (KOREAN 
HILLSIDE  BARRICADE,  TAEGU,  BLDG.  230,  STACK  CONFIGURATION  3) 


Target-Standing  Man:  Iigury  Criterion-Kinetic  Energy  Threshold:  58  ft-lb 
Monte  Carlo  Option:  Replications  60:  Seed:  17351:  Site  Altitude:  150  ft 
Stack  Top  Vertical  Fragmentation  Elevation  Zone:  70  to  110  deg 
90th  Percentile:  Stack  Barricade  (Hill)  Standoff:  32  ft:  Zero  Wind 


Range (ft) 


1/600 

Side  Fragnv  nation 
Top  Fragmentation 


FIGURE  14.  HAZARD  DENSITY  CONFIGURATIONS  FROM  SIDE  AND 
TOP  FRAGMENTATION  (KOREAN  HILLSIDE  BARRICADE,  TAEGU, 
BLDG.  230,  STACK  CONFIGURATION  3) 


28 


NAVSWCMP  89-353 


Hazard  Densities  for  Vertical  Fragmentation  Elevation  Zones  70-110 
and  80-90  deg  (Figure  15) 


In  the  fragmentation  arena  tests,  only  one  interaction  area  was  involved.  In  de¬ 
riving  the  top  fragmentation  from  the  side  fragmentation,  we  do  not  know  the  eleva¬ 
tion  zone  width  for  two  or  more  interaction  areas  making  up  the  top  fragmentation. 

The  two  elevation  zone  widths,  shown  in  Figure  15,  are  the  best  estimates  of  the 
maximum  and  minimum  zone  width.  From  Figure  15,  the  selection  of  the  70-  to  110- 
deg  zone  width  for  the  baseline  configuration  should  err  on  the  side  of  safety.  The 
reason  that  the  two  curves  come  together  at  approximately  1400-ft  range  is  tnat 
beyond  900-ft  range  only  the  side  fragmentation  is  involved  and  the  total  hazard 
density  beyond  900  ft  is  therefore  independent  of  the  top  fragmentation. 


Target  -Standing  Man  Injury  Criterion— Kinetic  Energy  Threshold:  58  ft-lb 
Monte  Carlo  Option:  Replications  60:  Seed:  17351:  Site  Altitude:  150  ft 
90th  Percentile:  Stack  Barricade  (Ilill)  Standoff:  32  ft:  Zero  Wind 


1/600 

70  -110*  (Baseline) 
80-90* 


*  Vertical  Fragmentation 
Elevation  Zone  (deg) 


FIGURE  15.  HAZARD  DENSITIES  FOR  VERTICAL  FRAGMENTATION  ELEVATION 
ZONES  OF  70-110  AND  80-90  deg  (KOREAN  HILLSIDE  BARRICADE, 
TAEGU,  BLDG.  230,  STACK  CONFIGURATION  3) 


Probability  of  Hazard  Hit 


The  hazard  density  threshold  (one  hazardous  fragment  per  600  ft2)  was  selected  to 
produce  a  threshold  probability  of  hazard  hit  of  approximately  0.01  when  the  pre¬ 
sented  area  of  the  personnel  target  was  6  ft2.  In  this  investigation  with  the  hill,  the 
angle  of  fall  for  all  fragments  is  very  steep  and  the  presented  area  of  the  standing 
man  is  much  less  than  6  ft2;  it  is  more  nearly  1  ft2.  With  the  hill  and  with  zero  wind,  a 
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rough  estimate  of  probability  of  hazard  hit  may  be  obtained  by  multiplying  hazard 
density  by  0.65.  If  a  probability  of  hazard  hit  criterion  were  used,  the  results  would 
have  appeared  much  less  hazardous.  If,  however,  we  had  used  a  prone  man  instead  of 
a  standing  man,  then  the  results  would  have  appeared  pretty  much  the  same. 


Positioning  the  Stacks 

In  terms  of  top  fragmentation,  there  is  very  little  effect  due  to  the  positions  of  the 
stacks  relative  to  the  hill.  In  terms  of  side  fragmentation,  however,  there  is  some 
effect.  Positioning  the  stack  near  the  front  of  the  shed  (BS  —  15  ft)  is  better  than 
positioning  it  at  the  back  of  the  shed  (BS  =  32  ft).  There  is  little  or  no  advantage  in 
positioning  the  stack  outside  the  shed  nearer  to  the  hill. 


CONCLUSIONS 


The  following  conclusions  may  be  drawn  from  the  results  of  this  investigation. 

The  pallets  of  Mk  82  bombs  should  be  stacked  in  a  compact  configuration  like 
Configuration  3,  shown  in  Figure  6. 

The  stack  should  be  positioned  as  near  to  the  front  of  the  shed  as  practical.  There 
is  very  little  advantage  in  positioning  the  stack  outside  the  shed  and  nearer  to  the 

hill. 

The  baseline  stack  Configuration  3  provides  a  significant  margin  of  safety  (see 
Figure  8)  under  the  hazard  criteria  currently  specified  by  DDESB.  At  800-ft  range, 
the  margin  of  safety  is  about  18  for  zero  wind. 

The  presence  of  tailwinds  do  not  significantly  alter  the  margins  of  safety  (see 
Figure  8).  The  tailwind  curves  are  still  well  below  the  threshold  hazard  of  one 
hazardous  fragment  per  600  ft2. 

Since  the  investigation  tended  to  err  on  the  side  of  safety  when  uncertainties 
arose,  hazards  are  not  expected  to  exceed  those  shown  in  Figures  8  through  15. 

The  hill  does  an  excellent  job  of  stopping  the  fragmentation  from  the  side  of  the 
stacks,  especially  when  the  stacks  are  positioned  at  the  front  of  the  shed.  The  hill  has 
little  or  no  effect  on  the  fragmentation  coming  from  the  top  of  the  stacks. 

In  the  Air  Force  request,  contained  in  the  appendix,  it  was  stated  that  the  storage 
shed  was  surrounded  on  three  sides  by  a  hill  with  the  fourth  side  open  and  blocked  by 
an  undefined  barricade.  This  fourth  side  should  be  investigated  when  detailed  data 
on  the  stack-barricade  geometry  become  available.  It  is  possible  that  a  compromise 
stack  configuration  may  be  necessary  to  minimize  the  overall  hazard. 
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